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ABSTRACT 



Aims. We explore the growth and the evolution of the bar instability in stellar-gaseous disks embedded in a suitable dark matter halo evolving 
in a fully consistent cosmological framework. The aim of this paper is to point out the impact of different gas fractions on the bar formation, 
inside disks of different disk-to-halo mass ratio, and the role of the cosmological framework. 

Methods. We perform cosmological simulations with the same disk-to-halo mass ratios as in a previous work where the gas was not taken into 
account. We compare results of the new simulations with the previous ones to investigate the effect of the gas by analysing the morphology of 
the stellar and gaseous components, the stellar bar strength and the behaviour of its pattern speed. 

Results. In our cosmological simulations, inside dark-matter dominated disks, a stellar bar, lasting 10 Gyr, is still living at z = even if the 
gaseous fraction exceeds half of the disk mass. However, in the most massive disks we find a threshold value (0.2) of the gas fraction able to 
destroy the bar. 

The stellar bar strength is enhanced by the gas and in the more massive disks higher gas fractions increase the bar pattern speed. 
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In a pioneering paper lBerentzen et al. I (fl998h presented the nu- 
merical results describing the impact of a gas fraction inside a 
stellar disk on the stellar orbits and on the whole bar instabil- 
ity. Their model galaxy containing stars and gas was compared 
with a pure stellar model, subject to the same initial conditions, 
chosen to produce the development of a strong stellar bar. 

They showed that dynamical instabilities become milder in 
the presence of the gas component, and therefore a weaken- 
ing of the bar itself is observed as a gradual process. Moreover, 
the overall evolution is accelerated because of the larger central 
mass concentration and the resulting decrease in the character- 
istic dynamical time. 

The gas behaviour in disk galaxies and its connec- 
tions with the bar feature was stud ied in several papers ( 

IBerentzen et all d200lh: 
( 2003}): iBournaud et al ' 



Friedli & Benz 
Athanassoula I 



(2005) , 

(2006) ). 



' (Il993j)- f 
(2002) 



Fuxl dl999b : 
Bottema I 



Stinson et al 



Mich el-Dansac & Wozniakl d2006l) 
However in all these works the evolution of the 
disk or of the disk+halo simulated system arises in a iso- 
lated framework, outside the cosmological scenario, even 
if progressive efforts for improving the model, accounting 
for more recent knowledges coming from the cosmological 
hierarchical clustering scenario of structures formation about 



density distribution and concentration of the dark matter ( DM) 
halos, h ave been perform e d in th e last years dCurir & Mazzei 



d 1999b : iMazzei & Curir I d200lb : lAmanassouTa& Misiriotis 



(2002)). Here, for the first time, we analyse the growth and the 
evolution of the bar instability in stellar-gaseous disks embed- 
ded in a DM halo evolving in a fully consistent cosmological 
scenario. Our model cannot be viewed as a general galaxy 
evolution model, since the gradual formation and growth of 
the stellar disk has not been taken into account. Thus our 
work canno t be c o mpared with the more r ecent papers by 



Abadi et all d2003b. iGovernato et ail d2004 . iRobertson etal 



(2004), ISommer-Larsen et all ( 2006b where the formation 



Send offprint requests to: A. Curir 



of a realistic disk galaxy within the hierarchical scenario of 
structure formation in A -dominated cosmologies has been 
followed self-consistently. However our approach allows us 
to vary parameters like the disk-to-halo mass ratio and the 
gas fraction inside the disk, and to analyse the growth of the 
bar instability and its dependence on such parameters in a 
self-consistent cosmological frame work for the first time. 
Following the pioneering study by IBerentzen et all d 19980 . we 
compare the result s of such new set of simulations with those 
of a previous set dCurir et al. , 2006) including a pure stellar 
disk in the same scenario and with the same initial conditions. 
The plan of the paper is the following: in Section 2, we 
summarise our recipe for the initial disk + halo system fully 
described in Paper I. In Section 3 we present the cosmological 
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simulations and in Section 4 we point out our results. Section 
5 is devoted to our discussion and conclusions. 



2. Method 

We embed a gaseous and stellar disk inside a cosmological halo 
selected in a suitable slice of Universe and follow its evolu- 
tion inside a cosmological framework: a ACDM model with 
Q„, =0.3, Q. A =0.7, cr 8 =0.9, h =0.7, where Q. m is the total 
matter of the Universe, Qa the cosmological constant, cr g the 
normalisation of the power spectrum, and h the value of the 
Hubble constant in units of lOO/i -1 kms^'Mpc -1 . 

A detailed description of our method to produce the cos- 
m ological scenario where the disk is evolved has been given 
in 



und erdense env i ronme nt. From its accretion history (see Fig. 
1 in ICurir et al.1 (120061) ") we conclude that our halo undergoes 
no significant merger during the time it hosts our disk, nor 
immediately before. The halo density pr ofile is well-fit t ed by a 
Navarro, Frenk and W hite (NFW) form dNavarro et all dl996l) : 
Navarro et al.1 (119971) ') at z < 2. The concentration, C vi >, equal 
to Rvir/Rs, takes an high value, 18.1, confi rming that this 



halo d oes "form" at quite high redshift ( e.g. IWechsler et al 



(2002) for a discussion about the link between concentration 
and assembly history of the halo). The dimensionless spin 
parameter of the halo is 0.0 4 at z=2, near to th e average one 



for our cosmological model (Mailer et al., 2002). 



Curir et al. ( 20061). Here we prese nt a short overview of our 2.2. The baryonic disk 



recipe. Moreover, Curir et al. ( 2006h showed that the numeri- 
cal resolution does not impair their main result: in pure stellar 
disks, long living bars are a 'natural' outcome of the cosmolog- 
ical scenario. 

2.1. The DM halo 

To select the DM halo, we perform a low-resolution simula- 
tion of a "concordance" ACDM cosmological model. The ini- 
tial redshift is 20. 

From this simulation we identify the DM halos at z=0 in 
the massQrange 0.5- S-IO 1 ^ -1 M , with a standard friends-of- 
friends algorithm. 

We select one suitable DM halo with a mass M~l0 ll h~ l M Q (at 
z =0). We resample it with the multi-mass technique described 
in lKlvpin et al. (2001 ). The particles of the DM halo, and those 
belonging to a sphere with a radius 4h Mpc, are followed 
to their Lagrangian position and re-sampled to an equivalent 
resolution of 1024 3 particles. The total number of DM particles 
in the high resolution region is 1216512 which corresponds to 
a DM mass resolution of 1.21 10 6 /z _1 M & . The high-resolution 
DM halo is followed to the redshift z=0. 

After selecting the halo and resampling at the higher 
resolution the corresponding Lagrangian region, we run the 
DM simulation, to extract the halo properties in absence of 
any embedded stellar disk. The mass of our halo at z=0, 
1.03 • I0 u h^ M Q , corresponds to a radius, R vir - 94.1k' 1 
Kpc, which entails 84720 halo particles. The nearest DM 
halo more massive than l0 w h~ l M Q is ~ 1900/i _1 Kpc away 
from the centre of our halo; the less massive one, having 
mass of 4.6 ■ 10 7 /; _l Mo, is ~ 215/2 -1 Kpc away. Moreover, 
the behaviour of the density contrast, 6, is monotonically 
decreasing with the radius, and 6 falls below the unity value 
at ~ 550/T 1 , 450/i~\ 350/i _1 physical Kpc away from the 
centre of our halo at z = 0, z = 1, and z = 2 respectively. 
Therefore, we conclude that the selected halo is living in an 



1 In the following, we refer to the mass as the virial mass i.e. 
that enclosed in a sp here with overdensity S = p/p cr j, = 178 • flj], 44 
dNavarro et aill200Cl) . 

2 Halos have been identified using the friends of friends algorithm 
with a linking length Z = 0.15, i.e. the mean interparticle distances, 
and with more than 8 particles 



The spatial distribution of particles follows the exponential 
surface density law: p = poexp-(r/Vo) where ro is the disk 
scale length, ro = Ah Kpc, and po is the surface central den- 
sity. The disk is truncated at five scale lengths with a radius: 
Rdisk = 20h~ l Kpc. To obtain each disk particle's position ac- 
cording to the a ssumed density d istribution, we used the re- 

1986b . We used 56000 star par- 



jection method dPress et al 



tides and 56000 gaseous particles to describe our disk. The 
(Plummer equivalent) softening length, the same for DM, gas, 
and star particles, is 0.5/U 1 Kpc in comoving coordinates. The 
stellar radial velocity dispersion, <x s , is assigned through a 
Toomre parameter Q. Q is initially constant at all disk radii 
and it is defined as Q - where k is the epicyclic fre- 

quency, and £ the surface density of the disk. We plug a value 
of Q: 1.5, which corresponds to a warm disk. The gas parti- 
cles share the same velocity field as the stars. The constraint 
Q g > 1 where Q g = '"^y * , needed for the gaseous disk stabil- 



ity dFriedli & Benz 1119931) is always verified. 



We embed the disk in the high resolution cosmological sim- 
ulation, at redshift 2, in a plane perpendicular to the angular 
momentum vector of the halo and in gravitational equilibrium 
with the potential. Its centre of mass corresponds to the the 
minimum potential well of the DM halo. The initial redshift 
corresponds to 10.24 Gyr down to z = in our chosen cosmol- 
ogy. 



3. Simulations 

We performed eight cosmological simulations of a disk+halo 
system. We also performed one simulation in an isolated frame- 
work, using a NFW halo, to disentangle the effects of the dy- 
namical state of the halo, and in particular of the evolution- 
ary framework, on the growth of the bar instability (§ 14. j} . 
W e exploited the p ublic parallel N-body treecode GADGET- 
2(Springell|200l. 



3 We note, however, that C NFW is defined against i^oo. the radius 
enclosing a sphere with overdensity equal to 200 times the critical 
density of the Universe, and not against R rjr as here; therefore, in our 
cosmological model, it is always C NFW < C w >. At z = our halo has 
Cnf w ~ 14. 
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Table 1. Simulations: initial values 



N M disk gas fraction M DM R DM ar m 1/2 halo 



cl 


0.33 


0.1 


0.64 


1.94 


1. 


1.08 


cosmological 


c2 


0.33 


0.2 


0.64 


1.94 


1. 


1.08 


cosmological 


c3 


0.33 


0.4 


0.64 


1.94 


1. 


1.08 


cosmological 


c4 


0.1 


0.1 


0.64 


6.4 


0.9 


1.68 


cosmological 


c5 


0.1 


0.2 


0.64 


6.4 


0.9 


1.68 


cosmological 


c6 


0.1 


0.4 


0.64 


6.4 


0.9 


1.68 


cosmological 


cl 


0.1 


0.5 


0.64 


6.4 


0.9 


1.68 


cosmological 


c8 


0.1 


0.6 


0.64 


6.4 


0.9 


1.68 


cosmological 


il 


0.33 


0.2 


0.95 


2.87 


0.85 


1.5 


NFW 



I col: simulation number and simulation type (c: cosmological simulation, i: isolated simulation) 

II col: mass of the disk in code units (i.e. 5.9 x 10 10 MO) 

III col: fraction of gas, i.e. gas-to-disk mass ratio 

IV col: DM mass inside the disk radius in code units 

V col: halo-to- disk mass ratio in side the disk radius 

VI and VII cols: Efst athiou et al J ( 119820 parameters: a = r ~\ v m is the maximum rotational velocity, r„, the corresponding radius 
VIII col: type of halo used 



Table 2. Simulations: final results 



N 


M OM 




6 


Q b 


"max 


bulge 


bars in bars 


cl 


0.77 


2.39 


0.68 


0.43 


8.4 


y 


n 


c2 


0.78 


2.4 


0.1 


0.01 


n 


n 


n 


c3 


0.78 


2.4 


0.07 





n 


n 


n 


c4 


0.73 


7.43 


0.58 


0.32 


5.8 


n 


y 


c5 


0.73 


7.44 


0.6 


0.3 


5.4 


n 


n 


c6 


0.73 


7.56 


0.46 


0.23 


5.6 


n 


n 


c7 


0.73 


7.6 


0.46 


0.2 


5.4 


n 


n 


c8 


0.73 


7.6 


0.42 


0.21 


5.8 


n 


n 


il 


1. 


3.03 


n 


n 


n 


n 


n 



I col: simulation number and simulation type 

II col: DM mass inside the disk radius in code units 

III col: halo-to-disk mass ratio inside the disk radius 

IV co l: maximum ellipticity at z= 0: strong bar (Mazz ei & Curir L 
1200 ll) require e > 0.4 

V col: bar strength according to lCombes & Sanders I J 1 98 lh : stronger 
bars correspond to higher Q b values 

VI col: major axis (physical Kpc) corresponding to the maximum bar 
strength 

VII col: morphology of the inner region of the disk 

VIII col: peculiar features inside the disk 



The simulations run on the SP4 and CLX computers lo- 
cated at the CINECA computing centre (BO, Italy; grant 
inato003) and on OATo Beowulf-class cluster of 32 Linux- 
based PC (16 processor AMD K7, clock 700 MHz, and 16 
Pentium 4, clock 1700 MHz) at the Osservatorio Astronomico 
di Torino. The main parameters and the initial properties of our 
set of simulations are listed in TableQ] 

A global stability criteria for th e stellar bar inst a bility in a 
disk galaxy is the one analysed in lEfstathiou et alj dl982l) . In 
such a paper the parameters a r m and 



(^Jttt (where v m is the 
maximum value of the disk rotational curve, r m the correspond- 
ing radius, a — rp" 1 an d M is the disk mass) have been defined. 
Efstathiou et al.l dl982l) stated the criterion , .I'",,,, > 1.1 over 



mation. The values of these parameters are reported in Table 
Q] Simulation il corresponds to the evolution of the same disk 
as in simulation c2 but we embedded it in a NFW halo with 
the same virial mas s and particle num ber as in our cosmologi- 
cal halo at z=0 (see lCurir et al. (2006) for more details). Such a 
halo is spherical, isotropic and in gravitational equili brium. The 



proced ure used to construct this halo is described bv lHernquist 
(Il993l) . 



All the results in Table Q] can be compared with those of 
the homologo us simulations per formed with a pure stellar disk, 
presented in Curir et al. (2006). We repeated, indeed, one of 
the simulation of lCurir et al.1 d2006l) with the GADGET-2 code, 
since the time-stepping criterion we used is not present any 
more in the current version of the code. A comparison with the 
previous run does not show any significative change. In such a 
paper we also performed several numerical tests on the stability 
of the results. Here we summarise our findings: 

- a simulation performed with a mass resolution increased 
by a factor 8, i.e. with a force resolution increased by a fac- 
tor 2 (for a disk-to-halo mass ratio 0.1; see section A.2 
in Curir et al. 2006), gives rise to a long living bar (10 
Gyr old) as in the case of the lower mass and force reso- 
lution; note that in this test the mass and force resolution 
are respectively ~ 2 • IOVMq and e * 250/T 1 

pc, thus comparable to the resoluti on of recent works on 
galaxy formation simulations, e.g. dSommer-Larsen et al 



2006tlStinsonet al., 20061 iGovernato et all 12006): 



(o-MG) 1 ' 



the range 0.1 < ar m < 1.3 for a disk being stable to bar for- 



- embedding adiabatically a classically stable stellar disk in 
a isolated DM halo having a NFW profile, while keeping 
the mass and force resolution constant, brings to no bar 
instability, thus ruling out the possibility that the embed- 
ding procedure is the main driver of the bar instability (see 
Section 5 in Curir et al. 2006); 

- varying the stellar softening, at fixed mass resolution, from 
0.36fr'Kpc to 0.65h _I Kpc does not result in any signifi- 
cant changes as far as the bar instability is concerned (see 
section A. 3 in Curir et al. 2006); 
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- the same holds when the stellar mass resolution is reduced 
by a factor of 1/6 and increased up to 10 times the value 
adopted in our set of simulations (see A.2. in Curir et al. 
2006). 

These results guarantee that, in the dynamical range currently 
reached with tree-code N-body simulations, no mayor numer- 
ical shortcoming is expected, at least under the point of view of 
the gravitational evolution. Moreover, as a final test, we re-run 
simulation c2 using two times more gaseous particles, 112000 
instead of 56000. We find that the final bar dissolution and the 
gas inflow are the same as in the previous case. 



4. Results 

In this section we present the final ( i.e. at z=0) isodensity con- 
tours for stars and gas of simulations in Table Q] All such fig- 
ures have been built up with the same box-size, number of lev- 
els, and density contrast (see caption of Fig. [T}. 
We define as a measure of the bar strength the value of the el- 
lipticity, e = 1 - b/a (Table 13; a strong bar corresponds to 
e > 0.4. 

A more dynamic al measure of the bar stren gth at radius R 
has been defined by ICombes & Sanders I dl98ll) by using the 



parameter: Q,(R) = J^g: where Fj ax (R) = [d®(R,0)/d6] max 
is the maximum amplitude of tangential force at a given radius 
R and < Fr(R) >- R{d<&o/dR) is the mean axisymmetric radial 
force derived from the m - component of the gravitational 
potential at the same radius. The maximum value of Q,{R) pro- 
vides a measure of the bar strength Qi, for the whole galaxy. 
We evaluated the components of the gravitational force directly 
from the masses and the radii of the system produced by the 
simulations. We then represented the matrix Q t (i, j) - <^(,''j)> 
on a suitabl e two dimensional grid in analogy with the method 



outlined bv lButa & Block I ( 120011) . Naming <2m the maximum 
value of Q,(i, j) in quadrant k of the map, according with 

Buta& Block! d200ll) . we define: Q b = 



We listed the final values of the parameter Qb in Table [2] 
The values of the bar strength derived with both methods are 
consistent. 

In all our cosmological simulations, with the exception of 
simulations c2 and c3, a stellar bar is still living at z=0. 

4.1. Star and gas morphologies 

The strongest and the longest bar arises from simulation cl 
(Fig. [TJ. In this case, the "geometrical" strength of such a bar, 
namely the ellipticity, 1 - b/a, evaluated through the isoden- 
sity plots, is higher (0.68) than the on e measured in the homol- 
ogous case (0.52, ICurir et aL I d2006h l. The strength measured 
through the Qb parameter is also higher than the corresponding 
one in the non-dissipative homologous case as shown in Fig. [2] 
which compares the maps of Q t (i, j) final values of the stellar 
bars. The enhancement of the stellar bar strength is due to the 
gaseous bar which is superimposed, i.e coupled, to the stellar 
one along all the evolution, (see Fig. [1] bottom panel). In the 
simulation c2, the gas develops a 'grand-design' spiral struc- 
ture, more evident than in simulation cl (Fig. 01. Such an ef- 



Fig. 2. Maps of the ratios of the axisymmetric radial and tan- 
gential forces in the disk plane at z=0 for simulation cl (top 
panel) and the analogous case without gas (bottom panel). Box 
size and spatial resolution are the same as in Fig.l. 

feet can contribute to destroy the st e llar ba r, in agreement with 
recent findings by iBournaud et alj 020051) . In the same simu- 
lation, c2, the gas produces a barred feature from the earlier 
stages, lasting until z ~ 1.2. After such a redshift, the gaseous 
central condensatio n contributes to the ba r destruction, accord- 
ing to the claim of lBerentzen et alJ dl998l) . These authors show 
that the growth of a gas mass concentration in the disk dis- 
solves the regular orbits in the stellar bar. In such a simulation, 
the stellar bar decreases its strength and disappears completely 
at z ~ 0.15, whereas in the c3 case, which corresponds to a 
larger initial gas fraction, the bar disappears at z ~ 0.6. Fig. [4] 
shows the gas and star morphologies of simulatio n c2 at z-0. 
Simulation c4 shows a stronger stellar bar than in ICurir et al 
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Fig. 1. Face-on, edge-on and side-on isodensity contours (from left to right) of simulation cl at z=0 (see text); top panel shows 
the stellar component, bottom panel the gaseous component. Spatial resolution is always 0.5/z~' physical Kpc and the box size 
is 40 times the spatial resolution. Contours are computed at 11 fixed levels ranging from 2 x 1CT 4 to 0.015 in term of density 
fraction of stars (gas) within the spatial resolution to the total star (gas) density in the map. 




Fig. 3. Gas isodensity contours, x-y projection, of simulation c 1 (left panel) and c2 (right panel) at z= 1 . 



d2006h ( ellipticity 0.58 from Table |2] instead of 0.3 from their 
Table 3) and bar-in bar features along the whole evolution until 
z=0 (see Fig. [5). Simulations from c5 to c8, which correspond 
to increasing gas fractions, show the same but less enhanced 
features during the system evolution, so that at z=0 no bar-in- 
bar appear, and for simulation c8 (Fig. |6]l the final bar is the 
weaker one. In this set of DM dominated disk simulations, the 
action of the increasing gas mass is not enough to destroy the 
bar. We note that, also in these cases, the gaseous bar is coupled 
with that of the stars from the beginning, however simulations 



with the larger gas fraction restore the gaseous disk earlier, so 
that gas and star distributions decouple from z=l (Fig.|7]i. 

In Fig. [8] we show the behaviour of the ellipticity and of 
the Qi, as a function of the increasing mass fraction in disks of 
different masses. It is remarkable the agreement we recover be- 
tween the trends: the dynamical evaluation of the bar strength 
and the geometrical evaluation of it through ellipticities. 

Fig. [9] compares the gas and star morphologies of the 
non-cosmological simulation, il, with those of the homolo- 
gous non-dissipative case (simulation i4 in Tables 2 and 3 of 
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Fig. 4. As in Fig.Q]for simulation c2. 




Fig. 8. Behaviour of the bar strength and ellipticity at z= for 
our set of cosmological simulations with increasing gas frac- 
tion. Qb (dotted line) and ellipticity (full line) of our more mas- 
sive disks (i.e. disk-to-halo mass ratio 0.33); Qb (dot-dashed 
line) and ellipticity (dashed line) of our less massive, DM dom- 
inated, disks (i.e. disk-to-halo mass ratio 0.1). 



Curir et al. (2006)) after 10.24 Gyr. In simulation il the bar 
forms approximately at the same time as in the pure stellar 
case. However the presence of the gas destroys the bar after 
7.5 Gyr instead of 8.3 Gyr as in simulation c2. Thus the cos- 
mological framework fuels a longer living bar, moreover, when 
the disk-to-halo mass ratio is 0.33, in such a framework we 



find a critical value of the gas fraction, 0.2, able to destroy the 
bar. 

In the case of DM dominated disks, all the simulations here 
performed show a bar feature still present at z=0. This is en- 
hanced for the lower gas percentage, 0.1, and slightly weakened 
for the high er ones, compared with the corresponding no dissi- 
pative case ( Curir et al. . 2006h . but the bar cannot be destroyed 
by the gas presence. 

4.2. Mass inflow and circular velocities 

The gas, initially distributed over all the disk surface with the 
same density distribution as the stars, evolves by increasing its 
concentration depending on the disk-to halo mass ratio and on 
its initial mass fraction. Both such parameters affect the life and 
the strength of the bar, as discussed above. 

One crucial parameter to monitorate the ability of the 
central gas condensation to destroy the bar, is indeed the ratio 
of the central gas mass over the total disk mass, that we show 
in Fig.fTOlandfTTI Fig.[l0] for DM dominated disk simulations, 
show that the rate of inflow at each radius decreases for higher 
gas fractions: from 65% to 37% at 2 Kpc and from 56% to 35% 
at 5 Kpc as we go from simulation c4 to c8. Fig. QT| for our 
more massive disks, show the same trend but with a stronger 
decrease: from 75% to 35% at 2 Kpc and from 66% to 41% 
at 5 Kpc as we go from simulation cl to c2. In particular, for 
simulation c2, where the gas fraction is able to destroy the bar, 
the accretion rate at the outer radii becomes quasi stationary 
after z — 1 (Fig. [TT] middle panel) since with the weakening 
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Fig. 5. As in Fig.[T]for simulation c4. 



of the bar, the gas inflow is no longer forced. In simulation 
c3, such a behaviour appears earlier (Fig.QT] middle panel) as 
well as the bar disappears earlier than in the c2 case. 
The circular velocity, v e ,, r , computed assuming a spheri- 
cal mass distribution and therefore having the expression 
(GM/r) l/2 where M is the mass of a component of the system, 
gives a measure of the absolute concentration of such a com- 
ponent. In Figs. [T2] and Qj] (left panels) we show its behaviour 
at z — for gas, stars and dark matter in simulations c4 and c6. 
In the same Figs, (right panels) we also show the behaviour of 
the v c irc = r{d<&/dr) evaluated from the force field of the stellar 
and DM components. The difference of these values from 
the previous ones is due to the deviation from the spherical 
isotropy of the mass distribution. In particular, the wide spread 
we observe for the DM component indicates a wide anisotropy, 
wider at large radii. 



The impact of a central mass concentra tion (CMC) on the 



bar evo lution has been discussed in detail bv lShen & Sellwood 



(2004). Their study is not directly comparable with ours, since 
the growth of the CMC is modelled with an analytic potential 
instead of being the result of gas inflow like in our simulations. 
In their work it was shown that the destruction power depends 
on the CMC, given as a fraction of the disk mass, and on the 
'softening' of its potential, r, t , a parameter which co ntrols the 



compactness of the CMC. IShen & Sellwood I (|2004) find that 
for a hard CMC (namely with r s very small) only few percent 
of the disk mass is able to destroy the bar, whereas for a soft 
CMC ( r s =0.3 Kpc) a fraction larger than 0. 1 plays the job. In 



our scenario, r s qualitatively translates into the radius, r ? , where 
we set the evaluation of the gas mass concentration. We deduce 
that for our more massive disks (i.e. when Rdm — 1 -94 in Table 
1) a gas mass greater than 9% of the disk mass inside r g - 2 
Kpc is needed to destroy the bar. However for DM dominated 
disks, (i.e. with = 6.4 in Table 1) a value of 18% of the 
disk mass in the same r g cannot destroy the bar. 

Therefore, our cosmological simulations point out that also 
the value of Rdm is a crucial parameter in this play. 

4.3. Bar pattern speed 

The bar pattern speed, £l p is the angular speed of the the bar- 
like density wave as viewed from an inertial frame. Here it is 
evaluated following the position angle of the major bar axis 
during the disk evolution, where t he major axis is selected with 



the sa me criterion El quoted in Mich el-Dansac & Wozniak 



(2006), namely as the radius where the ellipticity profile 
reaches a maximum. In Fig. [14] the values of Q p as function 
of the redshift are shown for the two sets of our simulations: 
the more massive disks and the lighter disks respectively. Since 
the behaviour of Q. p is rather noisy, we provided also the er- 
ror bars calculated with the running average method, i.e. they 
are obtained averaging three consecutive values and computing 
their standard deviations. This error bar is more representative 
of the fluctuations of the quantity Q p than of the real errors 
on the measurements. On the other hand, the uncertainty in the 
measure of the position angle of the semimajor axis is always 
less than 10°. Such a limit leads to a maximum error of 6% 




Fig. 7. Gas isodensity contours (top panel), x-y projection, of simulation c4 (left panel) and c8 (right panel) at z=l; bottom 
panels show that corresponding contours for stars. 
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Fig. 9. Isodensity contours of simulation il at 10.24 Gyr (to p and central pane ls) compared with the stellar isodensity contours 
of the same case without gas (bottom panel, simulation c4 of Curir et al. ( 20061) ') at z=0. 



in the evaluation of £l p . The error bar related to such a maxi- 
mum error is presented at the bottom of each panel in Fig. 14. 
For the more massive disks, the stronger decrease is observed 
in simulation cl. This is also the simulation with the longer 
and the stronger stellar bar (see Table [2j, usually coupled with 
the higher decrease o f the pattern speed during the evolution 



( Athanassoula , 2002). Simulation c2 shows the same trend, the 



decrease of Q p is slightly lower, and the bar disappears before 
z= 0, as in simulation c3. However in such a simulation, the 
pattern speed increases by decreasing the redshift. 

We derive that in the more massive disks the increasing gas 
fraction produces a speed up of the stellar bar pattern since the 
gas presence both shortens the bar and decreases its elliptic - 
ity. Moreover such a beh aviour is consistent with the findings 
of iBerentzen et al.1 dl998l) . i.e. an accelerated evolution of the 
whole system when the gas is included. 



On the other hand, the DM dominated disks show the 
strongest decrease of Q. p with the redshift. The lower dynami- 
cal support of these light disks compared to the more massive 
cases indeed enhances the role of the gas which slows down the 
pattern speed more quickly. 



5. Discussion and conclusions 

We presented eight cosmologi cal simulat i ons w ith the same 
disk-to-halo mass ratios as in Curir et al. (2006). In order to 
study the impact of the gaseous component, here we included 
and varied its percentage inside disks with different disk-to- 
halo mass ratios. 

Between the more massive disks, simulation cl, with the 
smaller gas percentage, shows a long lasting bar, 10 Gyr old, 
stronger than the one developed in the pure stellar case. In 
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Fig. 12. Left panel: circular velocities obtained assuming a spherical mass distribution of gas (red dashed line), stars (yellow 
continuous line) and DM (blue long dashed line) at z=0 for simulation c4. Right panel: circular velocities deduced from the radial 
force (upper curve: dark matter, lower curve: stars) at z=0 for simulation c4. 
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Fig. 13. The circular velocities as in Fig.[T2]for simulation c6. 



simulation c2, a gas percentage of 20% appears as a threshold 
value for the bar life as far as such disk mass is concerned. 
In such a simulation, indeed, the gas is able to destroy the 
bar at the very end of the evolution (z=0.15). By increasing 
the gas percentage, as in simulation c3, the bar disappears 
earlier than in c2 (at z=0.6). We find that a gaseous mass 
concentration equal to 9% of the total mass of the disk inside 
a radius of 2 Kpc, is a lower limit for the bar dissolution. In 
our more massive disks, where the baryons gravitational field 
is compelling with that of the DM halo (i.e. disk-to-halo mass 



ratio 0.33), we find a threshold value of the gas fraction, 0.2, 
able to destroy the bar. 

The gas shows a trend to be more widely distributed by 
increasing its fraction, for a given disk mass. Moreover, for 
the same gas fraction, the mass inflow is larger in more mas- 
sive disks. These effects contribute to weakening the bar which, 
however, in our cosmological framework, cannot be destroyed 
in DM dominated disks. In these cases, indeed, the cosmologi- 
cal evolution of the DM halo drives the growth of the bar insta- 
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Fig. 14. Evolution of the bar pattern speed with the redshift. Left panel: simulation cl (red full line), c2 (black dashed line), c3 
(magenta dotted line). Right panel: simulations c4 (blue full line), c5 (light blue dotted line), and c6 (green dashed line). The 
error bars on the curves are obtained with the running average method, and the error bar at the bottom of each panel represents 
the uncertainty in the measure of the position angle of the bar (see text). 



bility so that the presence of the gas appears as a second order 
effect on its evolution. 

The DM dominated disks show indeed a behaviou r whic h 
is strongly driven by the cosmology as in ICurir et al. (2006). 
In these cases : i) we do not find any value of the gas fraction, 
in the range 0.1-0.6, able to destroy the bar, ii) even a high 
value of central gas concentration does not succeed in dissolv- 
ing the bar. The bars in these disks are not a classical product 
of the self-gravity, which is very weak, or of angular momen- 
tum exchanges, since the disk rotates very slowly, but they are 
features that strongly depend on the dynamical state and evo- 
lution of the cosmological halo. Therefore the classical results 
emphasising the gas impact obtained outside the cosmological 
scenario ar e no lo nger applicable. This confirms the results of 
Curir et al. (2006), where in the DM dominated disks the bar 
feature is triggered and maintained by the cosmological char- 
acteristics of the halo, namely its triaxiality and its dynamical 
state. 

Therefore the whole set of cosmological simulations we 
present here is suggesting a play between three parameters 
which drive the bar formation and dissolution in cosmology: 
the gas fraction inside a suitable disk radius, r g , the r g value 
itself, and the halo-to-disk mass ratio inside the disk radius. 
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